pennState 


Department  of  Mechanical  &  Nuclear  Engineering 


High-Pressure  Combustion  Chamber  Dynamics 


Vigor  Yang 

The  Pennsylvania  State  University 
University  Park,  Pennsylvania  16802,  U.S.A. 
* Email:  vigor@psu.edu 


Presented  at 

International  Symposium  on  Energy  Conversion  Fundamentals 

Istanbul,  Turkey,  June  21-25,  2004 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

22  JUN  2004 

2.  REPORT  TYPE 

N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

High-Pressure  Combustion  Chamber  Dynamics 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

The  Pennsylvania  State  University  University  Park,  Pennsylvania  16802, 
U.S.A. 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADM001793,  International  Symposium  on  Energy  Conversion  Fundamentals  Held  in  Istanbul, 
Turkey  on  21-25  June  2005.,  The  original  document  contains  color  images. 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

uu 

18.  NUMBER 
OF  PAGES 

51 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


pennState 


Why  Supercritical  Combustion  Research? 
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•  most  booster  engines  operate  at  supercritical  conditions 

•  current  understanding  not  sufficient  to  support  design  optimization 
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Liquid  Rocket  Chamber  Conditions 
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Flow  Diagram  of  RD-170  Engine 
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-  Energia  booster  and  Zenit  first  stage,  up  to  10  flights. 

-  LOX/kerosene,  one  main  two  boost  turbopumps 

-  806  ton  thrust  (vacuum),  337  seconds  of  I  ,  O/F  ratio  of  2.63 

-  Chamber  pressure  250  bar,  turbine  inlet  pressure  519  bar  and  temperature  772  K 
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Worldwide  Efforts  on  Supercritical  Combustion  Research  (1/2) 


Department  of  Mechanical  &  Nuclear  Engineering  = 

•  Tamura  et  al.  /  NAL  (Japan) 

•  Mayer,  Oschwald,  Haidn,  etc.  /  DLR  (Germany) 

•  Habiballah,Vingert,  Grisch,  etc.  /  ONERA  (France) 
Candel  et  al.  /  Ecole  Central  Paris  (France) 

•  Woodward,  Pal,  Santoro,  etc.  /  Penn  State  (USA) 
Talley,  Chehroudi,  etc.  /  AFRL  (USA) 

Blevins,  Morris,  etc.  /  NASA  Marshall  (USA) 


Worldwide  Efforts  on  Supercritical  Combustion  Research  (2/2) 
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•  Oefelein  /  DoE  Sandia  Lab.  (USA) 

•  Bellan  /  NASA  JPL  (USA) 

•  Farmer  /  U.  of  Nevada  (USA) 

•  Habiballah,  et  al.  /  ONERA  (France) 

•  Yang  /  Penn  State  (USA) 
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Shadowgraph  Results  -  LN2  into  GN2 

Chehroudi  et  al.,  AIAA  99-0206 ,  AIAA  99-2489 
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pcr=  3.39  MPa,  Tcr=  126  K,  T^BOOK,  Tin  =  99~120K 
uin=  10-15  m/s,  D-  =  0.254  mm,  Re  =  25,000  ~  75,000 
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Characteristics  of  Supercritical  Fluid  Jet 


Mayer  et  al. 
AIAA  1996-2620 
TLN2  =  105  K 
r^GN2  —  300  K 

ULN2  —  m^S 

Din  =1.9  mm 


•  Thermodynamic  non-idealities  and  transport  anomalies 
in  transcritical  regime 

-  rapid  property  variations  -  large  density  gradient 

•  Diminishment  of  surface  tension  and  enthalpy  of  vaporization 

•  Pressure-dependent  solubility 

•  High  Reynolds  number 
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Effect  of  Pressure  on  Turbulence  Scales 
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•  Pressure  increases  from  1  to  102  atm,  Ret  increases  by  102 

•  Kolmogorov  microscale  r//l~  Ret~3/4  (decrease  by  1.5  order) 

•  Taylor  microscale  A/lt~  Ret~112  (decrease  by  1.0  order) 
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LES  Formulation  of  Supercritical  Fluid  Dynamics 
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•  Favre-filtered 
conservation 
equations 


•  Closure 
requirements 


•  Thermodynamic  and  transport  properties  z,  Cp ,  //,  Z,  Dim 

•  Subgrid-scale  turbulence  interaction  R ,  L ,  C 

•  Chemical  kinetics 
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Equations  of  State 
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Evaluation  of  Thermodynamic  Properties 


Department  of  Mechanical  &  Nuclear  Engineering 

•  Sensible  enthalpy:  h(p,T)  =  h°(T)  +  A hac ( p,T ) 

•  Internal  energy:  u(p,T)  =  u"(T)  +  Aul  JpJ) 

•  Specific  heat  Cp(p,T)  =  c;(T)  +  ACp  exc(p,T) 

A hexc,  Auexc ,  ACpexc  =  dense  fluid  corrections 
h°(T),  u°(T),  Cp(T),  =  values  in  dilute-gas  limit 


Pressure-explicit  type  of  EOS: 
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Viscosity,  pg/cm-s  Compressibility  Factor,  Z 
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Thermophysical  Properties  of  Nitrogen 
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compressibility  factor 


specific  heat 
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Droplet  Vaporization  and  Combustion 
in  Quiescent  and  Convective  Environments 
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•  Liquid  oxygen  (LOX)  droplet  vaporization  &  combustion  in  hydrogen 
and  water  5  <  p^  <  300  atm 

500  <Too<  2500  K 
50  <  D0<  300  pm 


•  Hydrocarbon  droplet  vaporization  &  combustion  in  air  and  oxygen 

5  <  poo  <  200  atm 
300  <T00<  2500  K 
100  <  D0<  1000  pm 


•  Unsymmetrical  dimethylhydrazine  (UDMH)  droplet  vaporization  and 
decomposition  combustion 

1  <  p^  <  180  atm 
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Spherical  Mode  (100  atm,  0.2  m/s;  t=610  ps) 


pennState  Breakup  Mode  (100  atm,  15  m/s;  t=170  ps) 
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Flow  and  Temperature  Fields 
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Px=  1 00  atm,  Tx=  1 000  K,  ux  =20  m/s,  T0=100  K,  d0=50  jam,  H/R=8 
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Effect  of  Pressure  and  Velocity  on  Droplet  Lifetime 
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Convective  Velocity,  m/s 

effect  of  ambient  pressure  on 

-  thermophysical  properties 

-  critical  mixing  state 

-  convective  heat  transfer 

effect  of  ambient  velocity  on 

-  convective  heat  transfer 


*  atmospherical  condition 

-  Ranz  and  Marshall’s  correlation 

TJ_  yr  ^Re=0  _ _ J _ 

Tft  Re=0  h  1  +  0.3  Re12  Pr1/3 

•  supercritical  condition 

-  LOX/hydrogen  system 

TJ_  ^Re=0  _ _ 1 _ 

r /;Re=0  h  1  +  0. 15634  Re1 1  Pr^0: 
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Supercritical  Fluid  Injection 
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pennState  Shadowgraph  Images  of  Cryogenic  Nitrogen  Injection 

Mayer  et  ah  AIAA  2001-3275 
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(poo=  6.0  MPa,  Tx=  300  K,  unl=  4.9  m/s,  Tnl=  132  K,  Dm=  2.2  mm) 
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Computational  Domain  and  Grids 
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•  Kolmogorov  microscale  ij/lt~  Re~3/4 

•  Taylor  microscale  l//~  Re~v 2 

•  3.4  <  pch  <10.0  MPa  and  Din  =  0.254  mm 

•  3<  ^<5  jam 


total  grids 

225x75  x72  =  1,215,000 

mean  grid  spacing  in 
near  injector  region 

A  =  5  pm 


pennState 


Density  Gradient  Field 


Department  of  Mechanical  &  Nuclear  Engineering  = 

(pM=  9.3  MPa,  1^=  300  K,  uin=  15  m/s,  Tin=  120  K,  Din=  254 


Vp|  (kg/m4) 


1.0000E+M  3.9401E+04  1.5524E+05  6.1167E+G5  2.4100E+06 


pennState 


Time  Evolution  of  Density  Gradient  Field 
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(P„=  6.9  MPa,  T .  =  300  K,  uin=  15  m/s,  Tin= 


120  K,  Din=  254  nm) 


pennState 


Snapshots  of  Density  and  Temperature  Gradient  Fields 


Department  of  Mechanical  &  Nuclear  Engineering 


(pOT=  9.3MPa,  1=  300K,  uin=  15m/s,  Tin=  120K,  t=  1.550ms,  Din=  254|um) 
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Most  Energy  Containing  POD  Modes  of  Axial  Velocity 
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(Poo=  9.3  MPa,  T  =  300  K,  uin=  15  m/s,  1^=  120  K,  Din=  254  pm)  | 
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Iso-Surfaces  of  Pressure  and  Density  Gradients 
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(p„=  9.3  MPa,  1  =  300  K,  uin=  15  m/s,  Tin=  120  K,  Din=  254  pm) 
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Effect  of  Pressure  ou  Density  and  Temperature  Fields 
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Tr  =  300K,  uin=  15m/s,  Tin=  120K,  Din=  254|um) 


p  (kg/m3) 


1.2 
0.6 
QS  0.0 
-0.6 
-1.2 


40  110  180  250  320  390  460  530 

p  =4.3MPa 


p  =  6.9MPa 


pennState 


Power  Spectral  Densities  of  Velocity  Fluctuations 
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(Poo=  9.3MPa,  T  =  300K,  uin=  15m/s,  Tin=  120K,  Din=  254pm) 


Frequency  (Hz) 
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Frequency  (Hz) 


Large  density-gradient  regions  act  like  a  solid  wall  that  amplifies  the 
axial  turbulent  fluctuation  but  damps  the  radial  one. 
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Vertex  Shedding  Frequency 
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(P,=  9.3MPa,  Tx=  300K,  uin=  15m/s,  Tin=  120K,  Dm=  254|am) 


Frequency  (kHz)  Frequency  (kHz) 
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Linear  Stability  Analysis  of  Real  Fluid  Jet 
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Annular  Flow  a 
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Density 
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Annular  Flow  a 
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Approach 

•  Two-dimensional  fluid  jet  instability 
at  supercritical  conditions. 

•  Unified  treatment  of  real-fluid 
thermodynamics  and  transport 
phenomena. 

•  Disperse  equation  solved  by  Newton- 
Ralpson  method. 


Conclusions 

•  As  the  density  ratio  increases,  the  spatial 
growth  rate  of  the  interfacial  instability  wave 
decreases.  Density  stratification  tends 
to  stabilize  the  mixing  layer. 

•  Density  stratification  has  little  effect  on  the 
frequency  of  the  most  unstable  mode. 
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Bi-Propellant  Swirl  Co-Axial  Injector 
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(ra)  j 
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Large  Eddy  Simulation  of  Swirling  Oxygen  Jet 

Department  of  Mechanical  &  Nuclear  Engineering  = 

Issues 


Swirling  jet  dynamics  at  supercritical 
conditions. 

Flame  stabilization  mechanisms  of  swirl 
co-axial  injector. 

Liquid  rocket  thrust  chamber  dynamics. 


x  (m) 


T(K) 


I 


Major  Results 

Liquid  film  thickness  and  swirl  cone  angle. 

Detailed  flow  structures,  including  central 
recirculation  zone,  surface  instability,  etc. 

Response  of  injector  dynamics  to  external 
forcing. 
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LOX/Kerosene  Preburner  Swirl  Injector 

(Wu,  et  al.,  unpublished  data,  2003) 
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Secondary  injection  Secondary  injection 


oxidizer-rich  preburner  injector 


damaged  inner  centrifugal  injector 


pennState 


Time  Evolution  of  Swirling  Jet 
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(p^lO.O  MPa,  T„=300  K,  uini=  30  m/s,  Tini=  120  K,  9=30°,  nitrogen  )  | 


Temperature  (K) 
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Time  Evolution  of  Swirling  Jet 
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(p^— 10.0  MPa,  1^=300  K,  Uj  =  30  m/s,  T-  •=  120  K,  0=30°,  nitrogen  ) 
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pennState  Disintegration  of  Swirling  Water  Jet 


(Inamura,  Tamura  and  Sakamoto,  JPP,  2003) 
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Water  Injection,  L/D=l  1.67,  K=1.0 


0  - 

30- 

mm 


(1)  M,  -26.1  g/s  (2)  M,  =32.5  g/s  (3)  M,  =48.6  g/s 

•  A  hollow  cone  sheet  forms  around  the  injector  exit. 

•  The  conical  sheet  fluctuates  vigorously  and  disintegrates  into  ligaments  and  droplets  at  the 
sheet  tip. 

•  The  sheet  breakup  point  approaches  the  injector  as  the  liquid  flow  rate  increases. 
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Theoretical  Analysis  of  Swirl  Injector  (1/2) 
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( Inamura  et  al. 
JPP,  2003  ) 
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Theoretical  Analysis  of  Swirl  Injector  (2/2) 
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Film  Thickness  at  Post  Exit 


C,  =  1 .429(1  +  A(x,  -  x' )}  - ().0279(x;  /  Re14 ) 


0  (deg.) 


Xt  <  X  <  XQ 


2.0 


1.5 


K  =  0.02798(x*  /Re1/4)  +  C: 


C,  =  1 . 1 43  -  0.02798(xl  /  Re1/4) 


£  i.o 


0.5 


x0  ={(1.1 82-  C2)/  0.2893}  Re 


1/4 


0.0 


L/de=1 1 .67 


•  • 


- it 


Mesurements 
-  Theory 


— U 


10  20  30  40  50 

M,  (g/s) 


pennState 


Limiting  Extremes:  2)  Diffusion  Processes  Dominate 
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"High"  Heating  Rates 

-  diminished  intermolecular  forces  promote  diffusion 
processes  prior  to  atomization 

-  injected  jet  vaporizes  forming  continous  fluid  in  presence 
of  exceedingly  large  interfacial  gradients 

-  diffusion  flame  resides  within  annular  post  wake 
separating  oxygen  jet  from  outer  hydrogen  flow 
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Burning  LOX  Jet  at  Supercritical  Pressure 

(Mayer,  DLR,  Germany ;  Tamura,  NAL,  Japan) 
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0*lox  ^  m/s,  Ujj^  300  m/s,  TLqX  100  K,  TH 


=  300  K,  p  =  6  MPa) 
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Supercritical  Injector  Flow  and  Flame  Dynamics 


Combined  OH  emission  and  backlighting  images  (Ph.D  thesis  of  Matthew  Juniper) 
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LOX/Hydrogen  Shear-Coaxial  Injector  Dynamics 
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Vorticity,  s'1  (x  10"6),  -1.5  <  co  <  0.4 
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Modeling  and  Simulation  of  Supercritical  Combustion 
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Effect  of  Pressure  on  Mean  Temperature  Distributions 
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<?„=  300  h  U|n=  15  ^  Tm=  120  K’  Din=  254  |im)| 


pennState  Effect  of  Pressure  on  Mean  Velocity  Distributions 
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(T„=  300  K,  uin=  15  m/s,  T,  =  120  K,  Din=  254  |im) 
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Frequency  Spectral  of  Radial  Velocity  Oscillations 
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p  =  4.2  MPa  I  P  =  6.9  MPa 


p  =  9.3  MPa 
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Normalized  Density  and  Temperature  Distributions 

along  Radial  Direction 


Department  of  Mechanical  &  Nuclear  Engineering  = 

T00=  300K,  uin=  15m/s,  Tin= 


120K,  Din=  254jum) 


•  Thermal  diffusivity  of  nitrogen  is  relatively  lower  in  the  region 
where  the  temperature  is  near  the  critical  temperature. 

•  Most  thermal  energy  transferred  from  the  hot  ambient  gaseous  nitrogen 
to  the  cold  jet  is  used  to  facilitate  volume  expansion. 
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Numerical  Challenges 
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Challenges 

•  machine  round-off  errors  at  low  speeds 

•  eigenvalue  disparity 

•  time  accuracy 

•  real-fluid  behavior 

•  robust  and  efficient  numerical  treatment 


High  Pres. 
Supercritical 
Gas  layers 


Solutions 

•  pressure  decomposition 

•  preconditioning  method 

•  dual  time-stepping  integration  technique 

•  partial  mass/molar  properties 

•  derivation  of  numerical  Jacobians  and 
thermodynamic  properties  based  on 
fundamental  thermodynamic  theories 
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Normalized  Density  and  Temperature  Distributions 
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(Pa=  9.3MPa,  T„=  300K,  uin=  15m/s,  Tin=  120K,  Din=  254pm) 


•  Due  to  the  “near  critical  slow  down”,  the  temperature  of  nitrogen 
fluid  increases  slowly  along  the  jet  centerline. 

•  A  self-similar  density  profile  exist  when  x/d>  15. 
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